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Abstract 

Evolution  of  vortices  and  their  interactions  are  investigated  in  the  Lagrangian  context  using  the  continuous 
vortex  sheet  method.  Among  the  topics  investigated  are  the  numerical  Helmholtz  instability  on  a  circular 
vortex  sheet,  the  fractal-pattern  instability  on  an  expanding  vortex  spiral,  interactions  of  vortices  with  shear 
layers,  and  some  other  phenomena  related  to  experimental  observations  in  the  laboratory  and  nature. 
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Subscripts 

n  normal 

t  tangential 

Superscripts 

+  to  the  left  of  C 

-  to  the  right  of  C 

*  complex  conjugation 


Fourier  coefficients 
oriented  contour 
Cauchy  density 

width  of  constant-density  segment 
imaginary  unit,  V-l 
positive  integer 
number  of  segments  of  C 
velocity 

velocity  components 

complex  velocity 

Cartesian  coordinates 

complex  coordinate  of  observation  point 

spanwise  circulation 

vortex  density 

Dirac  delta  function 

radius,  wing  half-span 

angle  between  real  axis  and  normal  to  C 

pseudo-time 

complex  coordinate  of  running  point 
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Introduction 

This  paper  studies  some  fundamentals,  although  in  a  simplified  form,  of  the  formation  and  evolution  of 
vortices.  The  stalling  point  is  a  two-dimensional  vortex  sheet  along  which  the  fluid  velocity  is 
discontinuous.  A  global  characteristic  of  the  vortex  sheet  evolution  is  stretching,  implying  that  the  vortex 
density  (strength)  of  the  sheet  is  a  decreasing  function  of  time.  Locally,  there  may  be  short-lived 
contractions  of  the  sheet  in  which  strength  increases  sharply,  preceding  sudden  qualitative  changes  such  as 
the  onset  of  roll-up  and  instabilities. 

The  majority  of  present-day  vortex  methods  have  evolved  from  the  point-vortex  approximation  by 
Rosenhead  [1],  A  schematic  of  such  a  scheme  is  shown  in  Fig.l.  The  vortex  particles,  which  are  created  at 
the  solid  boundary  at  each  time  step,  in  order  to  satisfy  the  no-slip  boundary  condition,  are  freed  in  the 
subsequent  time  steps  to  move  with  the  flow  and  diffuse  [2] .  However,  in  order  for  the  numerical  simulation 
to  be  accurate,  the  panicles  must,  to  a  certain  degree,  be  uniformly  distributed.  This  requires  the  usage  of 
techniques  such  as  redistribution,  splitting  and  merging  of  particles  [3].  Results  obtained  in  this  way  are 
equivalent  to  gridless  solutions  of  the  incompressible-flow  Navier-Stokes  equations. 

The  continuous  vortex  sheet  method  [4]  employed  in  this  paper  is  relatively  simple  in  comparison,  but  does 
not  have  the  capabilities  of  creating  or  diffusing  vorticity  (or  no  attempt  has  been  made  so  far).  It  does, 
however,  model  the  jump  discontinuity  across  of  thin  vortex  sheets,  Fig. 2,  preserving  discrete  spiral 
structures  of  vortices  similar  to  those  observed  in  flow  visualizations  [5] . 


Fig.  1  Discrete  vortex  representation 
(after  Leonard  [2]) 


Fig. 2  Continuous  vortex  sheet 

representation 


Continuous  sheet  model 


The  complex  velocity  induced  by  an  oriented  vortex  sheet  C  at  an  observation  point  z  is 

mo  i  f/(0 


w(z)  =  | 


2  TT(C-Z) 


\dO=  —  \ 

97 ri  J 


2ni  j,  £ 


-d£ 


(1) 


where  £  is  a  point  on  C  and  y  is  the  real-valued  vortex  density  function.  The  complex  contour  element  <f£ 
can  be  expressed  in  terms  of  the  arc  length  element  I  <f£  I  as 

d£=ld£lfeiv(C) 


where  v  is  the  angle  between  the  real  axis  and  the  normal  to  C,  pointing  to  the  right,  see  Fig. 2.  Accordingly, 
the  complex  density  function  in  Eq.(l)  is  given  by 

f(0=ir(Oe-iv(° 


(2) 
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The  second  integral  of  Eq.(l),  also  known  as  the  Cauchy-type  integral,  describes  a  function,  which  is 
analytic  everywhere  except  across  C.  In  accordance  with  the  Plemelj  formula, 

w+(0-w-(0  =  /(0  (3) 

where  the  superscripts  +  and  -  denote  the  limiting  values  of  w  to  the  left  and  right  of  C  respectively.  If  z  is  a 
point  on  C  ,  the  complex  velocity  is  calculated  by  taking  the  principal  value  of  the  Cauchy-type  integral  in 
Eq.(l).  Expressing  the  complex  velocity  on  C  in  terms  of  the  normal  and  tangential  components, 

w(C)  =  [?„(C)-^(C)]e“'V(0  (4) 

it  follows,  from  Eqs.(2)-(4), 

q+„(0-q~„(0  =  o 


and 

q;(0-q;(0=-no  (5) 

The  self-induced  time  evolution  of  C  is  calculated  by  integrating 

—  =  w\z)  (6) 

dr 

where  the  term  on  the  right-hand  side  is  the  complex  conjugate  of  the  principal  value  of  Eq.(l). 

The  integration  is  implemented  by  dividing  the  contour  into  n  straight-line  segments  (panels)  of  constant 
vortex  density,  and  advancing  the  solution  in  finite  steps  of  the  (pseudo-)time  T.  From  Eq.(l)  it  follows  that 
if  y  is  chosen  to  be  non-dimensional,  then  also  w  is  non-dimensional  and,  from  Eq.(6),  T  must  have  the 
same  dimension  as  z  (usually  normalized  by  span).  Conservation  of  vorticity  implies 

y(£)  I  dC,  I  =  constant 

Accordingly,  if  a  vortex  segment  stretches,  its  density  decreases  and  if  it  contracts,  its  density  increases.  A 
detailed  description  of  the  numerical  algorithm  is  given  in  Ref.  [4]. 


Vorticity  circuits 

As  an  example  we  consider  C  to  be  the  circle  of  radius  p  : 

Z=peiv 

A  variety  of  vortex  flows  can  be  generated  by  prescribing  the  vortex  density  as  the  Fourier  series 

7(v)  =  0.5<30  +  cos  jv  +  bj  sin  jv) 

j>  i 

where  y(v)  is  to  be  interpreted  as  y(v(£)). 


(7) 


Two  elementary  cases  will  be  analyzed  here:  1)  a  cylindrical  volume  of  fluid  moving  in  fluid  at  rest  and,  2)  a 
cylindrical  volume  of  fluid  at  rest  around  which  fluid  circulates. 

The  first  case  is  obtained  by  choosing 

y(v)  =  al  cosv  =  -y(c'1' +<?“'v)  (8) 
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2“ 
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From  Eqs.(2)  and  (7) 
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Substituting  in  Eq.(l)  and  using  the  Cauchy  formula  it  follows 


W(z)  =  i  y 
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j  r  -7  b  J  r  -v  b 
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This  result  can  be  interpreted  as  follows:  inside  the  circle,  there  is  a  uniform  downward  motion  described  by 
the  velocity  components 

u  —  0  and  v=-0.5flj 


Outside  the  circle,  the  velocity  is  that  of  a  point  doublet  of  strength  Kci\p 2 ,  oriented  in  the  vertical  direction. 
It  describes  the  velocity  of  ambient  fluid  displaced  by  the  moving  cylinder. 


The  evolution  of  the  vortex  circuit  whose  initial  conditions  are  given  by  Eqs.(7)-(8)  is  shown  in  Fig.3.  The 
left-hand  figure  column  shows  the  changes  of  the  circuit  geometry  and  the  right-hand  column  is  the 
corresponding  histogram  of  yin  the  50%  reduced  length  scale.  It  is  seen  that  the  circle  deforms  smoothly 
until  r  =0.72  (normalized  by  diameter).  At  that  instant  two  density  peaks  occur,  followed  by  the  onset  of  two 
lee-side  rollups.  The  sudden  appearance  of  singularities  (infinite-curvature  points)  in  the  shape  of  an 
evolving  vortex  sheet  is  a  phenomenon  well  explored  theoretically  [6]  [7].  The  subsequent  evolution  of  the 
vortex  circuit  in  Fig.3  is  as  expected:  a  continuous  stretching  driven  by  the  two  roll-ups. 


The  second  example  is  obtained  by  choosing 

y(v)  =  0.5a()  =  constant 


(9) 


From  Eqs.(2)  and  (7) 


no=iaff 


and  from  Eq.(l) 


w(z)  =  i 
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Inside  the  circle,  there  is  no  flow.  Outside,  the  flow  field  is  identical  with  that  generated  by  a  point  vortex  of 
strength  K  a0  p,  located  at  the  centre  of  the  circle.  Since 

q>  o 


it  follows  from  Eqs.(5)  and  (9) 

q~  =0.5a0 

The  circumference  of  the  circle  turns  with  the  average  of  the  above  velocities  (principal  value): 

qt  =0.5(q?  +  q~)  =  0.25a0 

As  may  be  conjectured  from  Fig.4,  the  circle  keeps  turning,  with  flow  swirling  around  twice  as  fast,  until 
about  T  =4.  At  this  (critical)  time  the  accumulated  numerical  errors  have  already  created  noticeable  wiggles 
in  the  shape  of  the  circle  and  small  variations  of  y  appear.  Subsequently,  a  rapid  breakdown  of  the  vortex 
sheet  takes  place:  multiple  peaks  of  ypop  up,  followed  by  a  series  of  roll-ups  on  the  circle  circumference, 
which  are  accompanied  by  an  intensive  stretching.  The  observed  long-term  existence  of  a  slightly  perturbed 
equilibrium  solution,  terminated  by  a  sudden  transition,  is  consistent  with  the  theoretical  analyses  [8]. 


The  difference  between  the  two  cases  shown  can  be  summarized  as  follows.  The  variation  of  yin  Fig.3 
allows  a  creation  of  density  peaks  early  in  the  vortex  evolution  process,  enabling  two  orderly  roll-ups  that 
facilitate  vortex  sheet  stretching.  In  Fig.4,  because  of  the  circular  shape  of  C  and  y  =  constant,  the  natural 
stretching  is  suppressed  -  until  the  Helmholtz  instability  takes  over. 
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Fig. 3  Evolution  of  vorticity  circuit, 

y=  cosv,  n  =  180;  At-  0.02 


Fig.4  Evolution  of  vorticity  circuit, 
y=0.5  ,n  =  180;  At  =  0.02 
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0.00 


Fig.5  Interaction  of  vorticity  circuit 
y=  cosv,  n  =180  and  shear  layer.  y=  -0.4,  in 
crosswind,  u  =  -0.2;  Az  =  0.02 


Fig. 6  Interaction  of  vorticity  circuit 
}fc0.5,  n=180  and  shear  layer,  y  =  -0.4, 
in  crosswind,  u  =  -0.2;  At=  0.02 


Interactions 

Vortex  sheet  roll-up  can  be  also  induced  externally,  for  example  by  an  interaction  of  several  vortex  sheets. 
In  the  first  frame  of  Fig.5  we  show  a  straight  vortex  sheet  of  constant  density.  It  is  combined  with  a  constant 
crosswind  in  such  a  way  that  there  is  no  flow  above  the  sheet  and  flow  to  the  left  below.  If  isolated,  the  sheet 
would  eventually  be  subject  to  a  Helhmholtz  instability  similar  to  that  developed  on  the  circle  in  Fig. 4. 
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However,  this  situation  changes  when  a  vorticity  circuit,  similar  to  that  shown  in  Fig. 3,  is  implanted  (at  time 
T  =  0).  The  circuit  moves  down  into  the  vortex  sheet,  first  indenting  it  and  then  inducing  vortices  at  the 
edges  of  the  indentation  basin.  The  interaction  is  not  symmetrical  because  one  of  the  impinging  circuit's 
vortices  is  of  the  same  sign  as  the  straight  vortex  sheet  whereas  the  other  is  of  opposite  sign.  This  flow 
pattern  resembles  the  interaction  of  the  counter-rotating  aircraft  vortices  with  an  atmospheric  shear  layer. 

The  example  in  Fig. 6  is  perhaps  even  more  interesting  because  it  combines  the  same  vortex  sheet,  with  no 
flow  above,  with  a  vorticity  circuit  of  Fig. 4,  with  no  translation  motion  of  its  own.  It  is  the  circulation 
around  the  vortex  circuit  that  initializes  the  interaction:  the  straight  vortex  sheet  is  pushed  down  on  the  left 
and  pulled  up  on  the  right.  The  deformed  vortex  sheet  no  longer  provides  zero  flow  above  and  thus  interacts 
with  the  vorticity  circuit.  Since  in  this  case  the  induced  shear  layer  vortex  on  the  right  rotates  flow  in  the 
direction  opposite  to  that  of  the  vortex  circuit,  both  systems  rise.  A  similar  effect  can  be  also  detected  in 
Fig.  5.  The  couple  of  equally  oriented  vortices  on  the  left  have  a  tendency  to  rotate  about  each  other  but  the 
couple  of  oppositely  oriented  vortices  on  the  right  translate  upwards. 

Vorticity  lines 

Besides  the  vorticity  circuits  it  is  also  of  interest  to  investigate  finite-length  stretches  of  vorticity.  Most 
familiar  are  the  trailing  vortex  sheets  of  lifting  wings,  known  to  roll  up  into  counter-rotating  wing  tip 
vortices.  In  essence,  the  integrals  of  Eq.(l)  are  singular  at  the  ends  of  an  open  contour  C,  providing  nuclei  of 
roll-up,  unless  the  vortex  density  is  vanishing  there.  Here  we  shall  consider  a  vortex  segment  parallel  to  the 
x-axis,  having  a  “bell-shaped”  density  distribution 

/(x)  = -^=exp(-k2x2) 

where  k  is  a  positive  integer.  It  can  be  shown  that 

|  y(x)dx  =  1 


and 


lim/(x)  =  S(x) 


where  5  is  the  Dirac  delta  function.  In  Fig.7,  we  immediately  notice  that  the  peaking  of  /during  the  vortex 
sheet  evolution  does  not  follow  the  limiting  process  k—>  °°.  Although  the  condition  of  constant  area  under 
the  /curve  is  satisfied  by  conservation  of  vorticity,  the  base  does  not  narrow  when  peaking  occurs,  leading 
to  a  “thorn-like”  /distribution.  In  any  case,  it  would  be  unrealistic  to  expect  that  a  cumulative  (singular) 
vorticity  could  be  generated  during  a  spontaneous  vortex  evolution  process. 


The  case  shown  in  Fig. 8  is  obtained  by  cutting  the  vortex  segment  in  half  and  inserting  a  segment  of  width 
h,  having  the  constant  density 


Y(x)  =  -=, 


\x\<h!2 


On  the  left  and  right  flanks  thus 


/(x)  =  —j=  exp[-k 2  (I  x  I  —h  /  2) 2  ] , 
dn 


(11) 


I  x  l>  /i  /  2 


(12) 
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Fig.7  Evolution  of  vortex  segment,  type  Sa, 
Eq.(10),  k-5,  n  =  240;  At=  0.002 


Fig. 8  Evolution  of  vortex  segment,  type  SBa, 
Eqs.(ll)  and  (12),  k  =  5,  h  -  0.2,  n  =  240; 
At=  0.001 


The  y distribution  described  by  Eqs.(ll)  and  (12)  is  still  continuously  differentiable,  but  the  impact  of  this 
modification  is  startling.  As  may  be  observed  in  Fig. 8,  the  peaking  takes  place  at  the  ends  of  the  constant- 
density  interval  and  later  spreads  to  several  locations  in  between.  Compared  to  the  smooth  roll-up  in  Fig.7, 
the  evolution  of  the  constant-density  segment  in  Fig. 8  is  “noisy”,  displaying  Helmholtz-type  instabilities. 
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Although  the  real-world  phenomena  have  more  than  two-dimensions  and  the  underlying  physics  is  much 
different,  similar  vortex  patterns  are  also  observed  on  a  grand  scale  as  hurricanes,  Fig.9,  or  galaxies,  Fig.  10. 
Edwin  Fiubble  classified  the  regular  galaxies  into  three  major  types:  spiral  (S),  barred  spiral  (SB),  and 
elliptical  (E),  as  shown  in  the  diagram  of  Fig.ll.  Spiral  galaxies  have  a  bulge  at  the  centre  and  a  disk 
containing  spiral  arms.  Barred  spiral  galaxies  are  spirals  that  have  a  bar  running  across  the  center  of  the 
galaxy.  Elliptical  galaxies  do  not  have  a  disk  or  arms  and  the  stars  within  move  irregularly  [9].  Just  how  the 
different  types  are  established  is  believed  to  be  a  matter  of  initial  conditions  of  the  amount  of  mass  and 
angular  momentum  in  the  condensing  gas  cloud,  and  of  interacting  with  nearby  galaxies  [10].  By  Fiubble ’s 
classification,  the  computer  simulation  in  Fig. 7  would  be  of  type  Sa  and  that  in  Fig. 8  of  type  SBa. 


Fig.9  Hurricane  Fran,  Sept.4,  1996, 
(NOAA  /  NASA) 


Fig. 10  Whirlpool  Galaxy  (M  51/NGC  5194), 
C.  Messier  (1773)  /  P.  Mechain  (1781) 
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Sb  Sc 

EO  E5 

-CO-KT* 

SBa 

SBb  SBc 

Fig.  1 1  Hubble’ s  Tuning  Fork  Diagram 


(SYA)  21-10 


Vortex  spirals 


The  stretching  of  the  vortex  sheet,  as  it  winds  into  a  spiral,  is  known  to  have  a  stabilizing  effect  [11],  [12]. 
Locally,  the  stability  of  this  convection  problem  can  be  characterized  by  the  displacement  of  the  vortex 
segment  relative  to  its  size  (Courant  number) 


I  AC  I 


When  using  a  constant  time  stepping,  initially  an  unstable  behaviour  may  be  observed.  Eventually,  the  spiral 
stabilizes  because  \AC,  I  increases  due  to  stretching. 


The  four  examples  in  Fig.  12  were  obtained  for  a  trailing  sheet  behind  an  elliptically- loaded  wing, 
represented  by  the  lifting  line  with  the  spanwise  circulation 

T(x)  =  (p2  -  x2)112 ,  \x\<  p 

where  p  is  the  half  span.  The  (non-dimensional)  vortex  density  distribution  is  obtained  as 

y(x)  =  -^-  =  (x/p)[l-(x/p)2rm 

ax 

Since  the  magnitude  of  y  grows  asymptotically  towards  the  wing  tips,  the  segment  endpoints  have  been 
distributed  as 

Xj  /  p  =  -  cos(;r  jin),  0  <  j  <n 


n=120 


n=480 


n=240 


n=960 


Fig.  12  Numerical  instability  of  the  near- field  vortex-sheet  spiral, 
At~  0.02,  100  steps 


All  the  spirals  shown  in  Fig.  12  were  obtained  after  100  steps  using  a  constant  pseudo-time  step  At  =  0.02 
(normalized  by  wing  span).  The  spiral  built  of  120  segments  appears  to  be  most  stable.  The  small-scale 
undulations,  developed  in  the  very  early  stage  of  the  roll-up,  smooth  out.  The  spirals  with  a  more  refined 
segmenting  develop  a  chaos-like  labyrinth  pattern  that  persists  in  the  later  stages  of  their  development.  This 
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phenomenon  is  also  known  as  the  fractalization  by  folding  [13].  Although  the  fractal  curve  of  Koch  is  often 
used  as  an  example,  the  folds  of  the  spiral  relate  more  closely  to  the  curve  of  Peano  [14].  The  latter  curve  is 
described  by  a  continuous  function  of  the  arc  length  that  fills  a  bounded  planar  domain  (rectangle)  and  in  the 
limit  of  infinite  refinement  has  nowhere  a  derivative  [15].  The  labyrinth  pattern  obtained  inside  the  spiral  of 
Fig.  13  describes  a  vortex  core  to  the  extent  the  formulation  by  Eq.(l)  permits:  the  tangential  velocity 
decreases  toward  the  core  centre  like  in  viscous  flow  and  it  is  difficult  to  distinguish  the  individual  winds  of 
the  spiral. 

It  is  worth  noting  that  undulations  of  the  spiral  arm  similar  to  those  shown  in  Fig.  13  have  been  observed 
experimentally  by  Pierce  [16],  in  airflow  around  plates  accelerated  from  rest,  see  Fig.  14.  Evidently,  there  is 
a  relationship  between  real  physics  and  the  computer  simulation.  As  may  be  expected,  the  connection 
between  the  artificial  viscosity  introduced  by  the  finite  time  stepping  and  true  viscosity  holds  the  clue. 


Fig.  13  Velocity  field  of  a  fractal  vortex  spiral 
(detail  of  Fig.  12,  case  n  =  960) 


Fig.  14  Vortex  sheet  roll-up  behind  a  plate  with 
a  convex  front  surface  and  flat  rear 
surface. 

Reprinted  from  the  article  by  Pierce  [16] 
with  the  permission  of  Cambridge 
University  Press. 


Concluding  remarks 

It  has  been  shown  that  a  number  of  physical  phenomena  related  to  the  evolution  of  vortices,  such  as 
stretching,  contraction,  instabilities  and  interactions,  can  be  modeled  in  a  simplified  form  by  the  continuous 
sheet  method.  The  numerical  examples  given  in  this  paper  were  produced  by  a  code  VSR  (for  Vortex  Sheet 
Roll-up),  written  in  Visual  C++  as  a  Windows  application. 
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